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Anomalous elastic softening of SmRu 4 Pi 2 under high pressure 
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The filled skutterudite compound SmRu4Pi2 undergoes a complex evolution from a paramagnetic 
metal (phase I) to a probable multipolar ordering insulator (phase II) at Tmi ~ 16.5 K, then to a 
magnetically ordered phase (phase III) at Tn ~ 14 K. Elastic properties under hydrostatic pressures 
were investigated to study the nature of the ordering phases. We found that distinct elastic softening 
above Tmi is induced by pressure, giving evidence of quadrupole degeneracy of the ground state in the 
crystalline electric field. It also suggests that quadrupole moment may be one of the order parameters 
below Tya under pressure. Strangely, the largest degree of softening is found in the transverse elastic 
constant Ct at around 0.5-0.6 GPa, presumably having relevancy to the competing and very different 
Griineisen parameters SI of Tya an d Tn. Interplay between the two phase transitions is also verified 
by the rapid increase of Tya under pressure with a considerably large Q of 9. Our results can be 
understood on the basis of the proposed octupole scenario for SmRu4Pi2. 

PACS numbers: 62.20.Dc; 71.30.+h; 72.55.+S 
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I. INTRODUCTION 



S1T1RU4P12 is one of the more interesting members of 
the family of filled skutterudite RT4X12 compounds (R 
= rare earths; T = Fe, Ru, Os; X = P, As, Sb) which 
shows a variety of physical properties including super- 
conductivity, metal-insulator (MI) transition, magnetic 
ordering and heavy fermions (HF) Q, 0, ■ The interest 
in SmRu4Pi2 is due to its metal-insulator (MI) transition 
at Tmi ~ 16.5 K, a subsequent antiferromagnetic (AFM) 
transition at Tn ~ 14 K 0] , and a strange H-T magnetic 
phase diagram 0, H| . Tn is obscure at lower magnetic 
fields; however, it is distinctly visible in several measure- 
ments like thermal expansion [(| even in a zero magnetic 
field. With increasing magnetic field, Tmi increases while 
Tn decreases, resembling CeBg, which has an antiferro- 
quadrupolar (AFQ) ordering and a subsequent magnetic 
ordering 0. For this reason, S1T1RU4P12 was initially 
thought to be an AFQ system; however, increasing doubt 
is being thrown on this view, stimulating interest in the 
order parameters (OPs). SmRu4Pi2 crystallizes in a cu- 
bic structure with space group Im3 (Tjj), like other filled 
skutterudite compounds. Magnetic measurements show 
that Sm is trivalent with total angular momentum J = 
5/2. Specific heat measurements suggest a crystalline 
electric field (CEF) scheme consisting of a Tg7 ground 
state quartet and an excited doublet at about 60 K 
in the Th symmetry 0, • This scheme is plausible for 
understanding other measurements including the elastic 
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constant 0, . Noticeably, the quartet with both 
orbital and magnetic degeneracy, is a key point for un- 
derstanding the various properties of SmRu4Pi2- Hydro- 
static pressures of about 1 GPa in this work are assumed 
to have no substantial effects on the CEF scheme. 

The AFQ scenario for SmRu4Pi2 is under question as 
to several aspects. First, almost no elastic softening is 
observed above Tmi |l| and the elastic anomaly at Tmi 
is not very large, unlike a typical AFQ ordering such as 
in CeB 6 [llj and DyB 2 C2 |12j. Second, by application of 
magnetic fields, the specific heat anomaly of AFQ order- 
ing in CeBg is enhanced, while that at Tmi in S111RU4P12 
changes only slighy; meanwhile, the anomaly at Tn in 
CeBfi is weakened, whereas it is enhanced in SmRu4Pi2 
0, El • Recently, the possibility of octupole ordering was 
proposed for the MI transition ^4[. This proposition is 
based on the appearance of elastic softening in phase II 
toward Tn and the indistinctness of the subsequent mag- 
netic ordering at Tn, because these facts suggest a new 
coupling of $Me ($f: probable octupole OP in phase 
II; M: dipole OP in phase III; e: elastic strain induced 
by ultrasound) and indicate breakdown of the time re- 
versal symmetry (TRS) in phase II. Breakdown of TRS 
is confirmed by muon spin relaxation (/xSR) |15| , nuclear 
magnetic resonance (NMR) 16] and Sm nuclear resonant 
scattering [l^ measurements. An AFQ ordering is non- 
magnetic and holds TRS, so is ruled out as the primary 
OP in phase II. On the other hand, an octupole order- 
ing is magnetic and breaks TRS, so is possible in phase 
II. Microscopic measurements also suggest appearance of 
magnetic dipole moment in phase II even in a zero mag- 
netic field, in addition to an unknown multipole ordering 
[Rj|l7j |. This is very different to the strong candidate for 
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octupole ordering material, NpC>2 [l8j, which has only 
one phase transition. 

An iso-structural com pou nd, PrRu4Pi2, also shows MI 
transition (Tmi — 63 K) ■ Its MI transition is accom- 
panied by no magnetic anomaly, whereas a distinct mag- 
netic anomaly is observed in SmRu4Pi2. Moreover, a 
structural distortion is observed in PrRu4Pi2 [13] but is 
absent in SmRu4Pi2- These facts indicate a magnetic ori- 
gin for the MI transition in SmRu4Pi2, unlike PrRu4Pi2. 
In fact, a charge density wave (CDW) due to Fermi sur- 
face nesting is plausible for interpreting the MI transition 
in PrRu4Pi2 |2lj| . Recent X-ray and polarized neutron 
diffraction experiments have shown that the Pr-ion site 
splits into two nonequivalent crystallographic sites and 
suggest two different CEF schemes below Tmi i22| , where 
the importance of the p-f hybridization in the MI tran- 
sition is also emphasized. 

Ultrasonic measurement under hydrostatic pressures 
is rarely performed, partly due to its difficulty We suc- 
ceeded, however, in obtaining stable ultrasonic echoes by 
selecting a suitable sample of sufficient length. By per- 
forming such measurements, we intend to gain an insight 
into the following peculiarities in SmRu4Pi2- (1) The T 6 7 
quartet ground state degenerates with respect to and 
Ts-type quadrupole moments; however, no distinct elas- 
tic softening is observed at ambient pressure. Application 
of hydrostatic pressures is expected to induce stronger 
quadrupole-strain and inter-site quadrupole interaction, 
which then affects elastic properties. Quadrupolar order- 
ing, if visible under pressure, is helpful for understanding 
the possible OPs below Tmi and the peculiar H-T phase 
diagram. (2) We intend to build a pressure-temperature 
(P-T) phase diagram from which to gain knowledge of 
the OPs of phase II and III. (3) One more point of inter- 
est is the very different Gruneisen parameters fl of Tmi 
and Tn, the latter being much larger than the former 
|23j. Widely different Gruneisen parameters suggest dif- 
ferent pressure dependencies of the two phase transitions, 
and prompted us to explore the physical properties under 
hydrostatic pressures related to this phenomenon. 



II. EXPERIMENTS 

A cylindrical polycrystalline SmRu4Pi2 sample 4.1 mm 
in length and 3.4 mm in diameter was employed. It is the 
same one as used for ultrasonic measurements in mag- 
netic fields and was prepared at high temperatures and 
pressures using a wedge-type cubic anvil high-pressure 
apparatus. Actually, a single crystal is preferred here, 
but the available single crystals currently are very small, 
with dimensions of only lxlxl mm, preventing us from 
obtaining stable ultrasonic echoes under high pressures. 

The relative change of sound velocity v was measured 
by means of a phase-comparison technique. This tech- 
nique measures the relative Av/v by detecting the rela- 
tive frequency change A/// while keeping the phase shift 
at zero and neglecting any changes in sample length Al / 1 



(to be explained later). Elastic constant C is related to 
sound velocity v by the equation C — pv 2 , where the den- 
sity p of SmRu4Pi2 is 5.921g-cm~ 3 . LiNb03 piezoelectric 
plate with fundamental resonance frequency of 5-20 MHz 
was used as the ultrasonic transducer to both generate 
and detect the ultrasound. The only experimentally ac- 
cessible ultrasound modes for a polycrystalline sample, 
i.e., longitudinal Cl propagating parallel to polarization 
and transverse Ct propagating perpendicular to polar- 
ization, were measured as a function of temperature un- 
der different hydrostatic pressures. The absolute values 
of Cl, Ct and bulk modulus Cb of SmRu4Pi2 at 100 K 
are 153, 24.6 and 120 GPa, respectively @. 

Hydrostatic pressure up to 1.15 GPa was generated 
using a Cu-Be piston-cylinder pressure cell and a Teflon 
capsule inside it. The cell has two layers, consisting of 
an outer cylinder and an inner one placed inside. The 
inner diameter and the length of the inner cylinder are 
8 mm and 40 mm, respectively. The sample with two ul- 
trasonic transducers bonded to it, together with a tin 
manometer, is placed in the Teflon capsule filled with 
pressure medium (a 1:1 mixture of Fluorinert FC 70 and 
FC 77). The measurements were carried out down to the 
temperature of liquid helium, 4.2 K, or down to *~1.5 K 
by pumping. The effective hydrostatic pressure inside the 
Teflon capsule was calibrated by measuring the supercon- 
ducting transition temperature of the tin manometer. 



III. EXPERIMENTAL RESULTS 
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FIG. 1: Temperature dependence of the longitudinal elastic 
constant Cl under various pressures. The solid lines are cal- 
culated curves by Eq. 0. The arrows indicate T so ft below 
which elastic softening appears. 

Fig. n shows the temperature dependence of the lon- 
gitudinal mode ACl/Cl under various pressures. The 
solid lines are fitted curves that will be explained later. 



3 



The Cl at ambient pressure increases on cooling down 
to Tmi without clastic softening. Application of hydro- 
static pressures induces a clear softening above Tmi, but 
its pressure dependence is very weak. This may be par- 
tially due to the relative hardness of Cl that includes 
bulk modulus Cb- The temperature T S oft, below which 
softening appears, shows a slight increase with increasing 
pressure. 
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FIG. 2: Temperature dependence of the transverse elastic 
constant Ct under various pressures. The solid lines and 
arrows denote the same as in Fig. Q 

Fig-Hlshows the temperature dependence of the trans- 
verse mode ACt/Ct at various pressures, together with 
the fitted curves. Ambient pressure evidences a weak 
elastic softening above Tmi up to 35 K. With increasing 
pressure up to 0.6 GPa, elastic softening becomes increas- 
ingly marked. Furthermore, T S oft increases from 35 K for 
GPa to about 80 K for 0.6 GPa. Strangely, upon further 
increasing the pressure by 0.7 GPa, the elastic softening 
fades again, while at the same time, T so f t also decreases. 
This strange pressure dependence presumably reflects the 
complexity of the phase diagram with two competing 
phase transitions. A kink structure is observed at around 
90 K in both Cr and Cl at ambient pressure, reminiscent 
of the off-center motion of the R ions in the pnictogen 
cage, namely "rattling" as in PrOs4Sbi2 0|. However, 
here it seems not due to an intrinsic origin because it is 
absent in single crystals. The kink becomes very weak or 
undetectable as soon as even a relatively weak pressure 
is applied; this is also not expected for a rattling motion. 

For both Cl and Ct, a distinct slope change in the 
background C°, is noticeable i.e., the lattice contribu- 
tion to elastic constant, before and after applying pres- 
sure. As for the reasons, the influence of pressure- induced 
sample length variation Al/l, which alters the relative ul- 
trasound velocity as follows, is possible. 

Av Af Aip Al 

v f <p I 



Eq. shows that in the phase-comparison technique, 
relative sound velocity Av /v is composed of three parts, 
e.g., the relative frequency A///, relative signal phase 
Aip/ip and sample length Al/l. A(p/(p is kept at zero by 
a feedback loop that compares the phase shift between 
the ultrasonic signal transported by the sample and the 
reference signal. Al/l is usually negligibly small. There- 
fore, this equation also explains why ultrasonic velocity 
v is measured as a shift in frequency. If Al/l induced by 
pressure accounts for the change in C° , the linear thermal 
expansion coefficient a(T) is anticipated to be a change 
over 5 x 10~ 5 K _1 under a pressure of about 0.3 GPa, 
in comparison with that at ambient pressure. This is 
at least an order of magnitude greater than the a(T) at 
ambient pressure for SmRu4Pi2 [g, making it impossible 
for Al/l to account for the change in C°. Although this 
peculiarity is not well understood, it does not affect our 
analysis of the elastic softening at low temperatures. 

The temperature where a steep elastic drop occurs is 
known as the MI transition temperature Tmi at ambient 
pressure. Similarly, all such anomalies observed under 
pressure are ascribed to MI transition. Although sig- 
nificant elastic softening is induced by hydrostatic pres- 
sure, the shape of the abrupt anomaly at Tmi shows no 
substantial change. In fact, electrical resistance under 
high pressures remains characteristic of MI transition up 
to at least 1.2 GPa, while metallic behavior appears at 
P> 3.5 GPa [25| . Tmi increases on applying pressure, 
in agreement with the observations in electrical resis- 
tance. This issue will be discussed later in terms of 
Griineisen parameter £1. Unfortunately, we can not de- 
tect an anomaly at Tn under pressure, as was successfully 
done in magnetic fields [jj- 

IV. DISCUSSION 

A. Pressure-induced elastic softening 

Symmetrized strains induced by ultrasound perturb 
CEF state via quadrupole-strain interaction grOrer- 
Taking the inter-site quadrupole interaction g' r Or{Or) 
into account, elastic constant Cr can be formulated as 
follows by the second derivative of Landau free energy. 

C r (T) = C«{T)- N9 ^ T l (2) 

Here C{!(T) is background elastic constant, N the 
number of R ions per unit volume, Xr(T) the corre- 
sponding strain susceptibility calculated using the CEF 
scheme, and gr and g' r are the coupling constants of 
the quadrupole-strain and quadrupolc-quadrupole, re- 
spectively. Cr(T) measures the diagonal (Curie term) 
and non-diagonal (Van-Vleck term) quadrupolar matrix 
elements, analogous to the magnetic susceptibility that 
measures magnetic dipole matrix elements |26j . Analyz- 
ing the experimental results by using xr{T), one experi- 
mentally obtains \gr\ and g' r . 
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At ambient pressure, no elastic softening was observed 
above T MI for Cl and only very slight softening for Op- 
Generally, this indicates a very weak quadrupole-strain 
coupling <?r for the r@7 quartet ground state. This is be- 
lieved to be a signature of complex multipolar ordering 
in phase II. On the other hand, a high degree of soften- 
ing is observed under hydrostatic pressures, especially for 
Ct, giving evidence of the degeneracy of the CEF ground 
state with respect to quadrupole moment. Hydrostatic 
pressure induces elastic strain £b — £xx + e yy + e zz that 
corresponds to bulk modulus Cb. £_b has no direct cou- 
ple with symmetry quadrupole moment Or, thus is not 
a direct reason for the enhanced elastic softening. 

A crucial point in the physics of the enhanced elastic 
softening is which of the two types of inherent quadrupole 
moments (I^ or T§) is responsible. The answer to this 
question may give direct insight into the OP symmetry 
of phase II; unfortunately, it can not be clarified based 
on the present results on a polycrystalline sample. In 
contrast to the symmetrized transverse modes in a sin- 
gle crystal like {C\\ — C\i)/2 and C44 that respectively 
corresponds to T3 and IVtype quadrupole moment, Or 
and Cl in a polycrystalline sample reflect an average of 
different elastic modes. The former consists of C44 and 
(Cn — Ci2)/2; the latter all the available modes including 
Cn and bulk modulus Cb- For the same reason, Eq. @ 
cannot be directly used to analyze elastic softening here. 
Nevertheless, as long as the Xr(T) is dominated by the 
Curie term as predicted in the present system, Eq. (0) 
can be reformulated as follows. 



C°=-0.000087*T+1 .058 



C r (T) 



(3) 



According to Eq. , the quadrupole-strain coupling and 
inter-site quadrupole interaction can be evaluated as an 
average. Furthermore, transforming Eq. to Eq. Q is 
helpful in actual analysis. 



T-T a 



C°(T)-C T (T) T°-T Q 



(4) 



Here Tq = t?p | { | Or I ) 1 2 indicates the inter-site 
quadrupole interaction. The difference of the two charac- 
teristic temperatures, -Tq = Npr\(^Pi\Ov\(Pi)\ 2 /C-S, is 
usually termed the J aim- Teller coupling energy £'tt|27|. 
and is a function of quadrupole-strain coupling. 

Before analyzing elastic softening, elastic background 
C°(T) should be carefully estimated. This process may 
affect the fitted values of Tq and Ejt to some degree. 
Here we employ the commonly used T-linear C . Fig. [3] 
presents the measurement data of Ct and two estimated 
C°(T) lines. All the curves at non-ambient pressures 
are superposed on each other and a common C°(T) is 
estimated for them. Significantly, we found the pressure 
dependences of Tq and Ejt to be almost independent of 
C°(T) as long as a common background is employed as 
done in this Figure. In fact, a shared slope of all Or(T) 
curves under non-ambient pressures at temperatures over 
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FIG. 3: Transverse elastic constant Ct(T) under various pres- 
sures that is superposed on each other. The two solid lines 
indicate estimated elastic background, described as C° = 
-0.00021 x T + 1.042 and C° = -0.000087 x T + 1.058, for 
GPa and other non-ambient pressures, respectively. 



100 K can be observed, supporting the common C°(T). 
This procedure can maximally reduce extrinsic influence 
in estimating the two characteristic temperatures, and is 
crucial for the later discussions. 
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FIG. 4: (C°-Ct) -1 as a function of temperature under various 
pressures. The solid lines are Curie term fits based on Eq. , 
which give values of Tq and E JT = T<2 — Tq. 



Fig.0]shows (C°-Ct) _1 as a function of temperature. 
The T-linear increase above Tmi conforms to the Curie 
term described by Eq. |0J, due to the T 3 and/or r 5 - 
type quadrupole moment of the ground state. The 
horizontal intercepts and the slopes of the linear fits give 
the quadrupolar interaction Tq and Jahn- Teller coupling 
energy £gr ■ For longitudinal Cl we performed the same 
procedure (not shown here). Note that the calculated 
curves shown in Fig. ^ [21 are based on the obtained Tq 
and Ejt values, and Eq. ©. 

Fig. [S] presents pressure dependence of the quadrupo- 
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FIG. 5: Quadmpolar interaction Tq as a function of pressure. 
Open circles denote data estimated from Cl and solid circles 
from Ct- 



100 




FIG. 6: Pressure dependences of (a) the temperature T ao ft 
below which elastic softening appears, and (b) the Jahn- Teller 
energy _Ejt = Tq-Tq. Open circles denote data estimated 
from Cl and solid circles from Ct- 



lar interaction, — Tq vs. P. Tqs are all negative in indi- 
cation of antifcrro-quadrupolar interaction. Tq is once 
enhanced by pressure; however, it is weakened again 
at around 0.7 GPa where elastic softening is weakened. 
Shown in Fig. are T so f t and Jahn- Teller energy Ejt 
as a function of pressure. Both of them depend on the 
same quadrupole-strain coupling and reflect the lattice 
instability. As expected, T so f t and Ejt show a similar 
pressure dependence. However, the two energies exhibit 
different pressure dependences for Ct and Cl; enhanced 
values around 0.5-0.6 GPa are found in transverse Ct, 
while no distinct change is observed in longitudinal Cl- 
Here we first discuss the pressure dependence of 
Tq and Ejt- Because Tq — g^\(ipi\Or\fi)\ 2 and 



Ejt = Ng^\((fi\Or\<Pi)\ 2 /C^, two possible explanations 
of their pressure dependences are available. (1) The 
two energies change with pressure due to changes in 
quadrupole moment K¥>i|0r|v»)l- (2) Their changes are 
due to g' T and gr, respectively. For case (1), one expects 
a similar pressure dependence of Tq and Ejt • This is dif- 
ferent to our observations. Moreover, a significant change 
of the CEF scheme and the consequent quadrupole mo- 
ment are not anticipated under the present hydrostatic 
pressures. The quadrupole Kondo effect may have an in- 
fluence on K^jlOrlvi)! but one would expect a screened 
quadrupole moment by applying pressure, not an en- 
hanced one. Thus, it seems likely that changes in g' T 
and gr can be ascribed to the characteristic energies un- 
der pressure. Because 3p~Tq is negative and g 2 ~EjT is 
positive as factors in Eq. 10, the rapid decline of |Tq| 
and enhancement of Ejt around 0.5-0.6 GPa will en- 
hance elastic softening, consistent with the experimental 
facts in Ct- The absence of enhanced |Tq| and Ejt in 
Cl within the present experimental accuracy remains an 
open question. Besides the hardness of Cl, which pre- 
vents a precise analysis of elastic softening, the difference 
between Ct and Cl may also suggest that the softening 
is due to one of the possible quadrupole modes, T3 or 
but not both. 



B. Increase of Tmi with large Griineisen parameter 



Our results show that Tmi increases with pressure 
(Fig. , consistent with observations of electrical resis- 
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tance [25J . It is also in agreement with thermal expansion 
measurements where lattice shrinkage is observed below 
Tmi El- The monotonic increase of Tmi with pressure 
contrasts strikingly with the unusual changes in Tq and 
Ejt, suggesting again that the MI transition is not a 
quadrupolar ordering. Comparatively, no obvious pres- 
sure dependence [2^ or a very slight increase [2^ was ob- 
served for the Tafq in the typical AFQ compound CeB 6 . 
Formula 10 defines the Griineisen parameter as volume 
derivative of the phase transition temperature. 



lnT, 



MI 



c B m 



MI 



lnV Tmi dP 



(5) 



Based on this formula, the Griineisen parameter for 
Tmi is estimated to be 9 (solid line in Fig.0. The pres- 
sure dependence of Tmi reported in ref. [2{| leads to an 
fl of 7, roughly in agreement with ours. It is significant 
to make a comparison between the r2(TMi) of S111RU4P12 
and PrRu4Pi2- Performing the same calculation using 
the reported bulk modulus Cr = 207 GPa [3(j and the 
pressure dependence of Tmi j3l| . we obtain f2(TMi) — 2 
for PrRii4Pi2. The JI(Tmi) for PrRu4Pi2 is a normal 
value while that for SmRu4Pi2 is large, again reflecting 
the essential difference between the two MI transitions. 

The considerably large f2 is reminiscent of the two very 
different Griineisen parameters estimated for Tmi and Tn 
[2^. The conventional equation Q — CbP/C (here /3 is 
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FIG. 7: Tmi as a function of pressure. The solid and open 
circles denote Tmi observed in Ct and Cl, respectively. The 
solid line is a guide for eyes and represents a Griineisen pa- 
rameter of £1 = 9. The broken lines represent two supposed 
pressure dependences for Tn (£1 = 28) and Tmi (£1 = 1), re- 
spectively. 



volume thermal expansion coefficient and C specific heat) 
gives a huge = 28 for Tn, but a normal £1 = 1 for Tmi- 
Widely different Os are expected to produce widely dif- 
ferent pressure dependences of T MI and T N , as indicated 
by the broken lines in Fig.0 Unfortunately, no magnetic 
fields were applied in this work, so Tmi and ?n are n ot 
separated and a more detailed P-T phase diagram is not 
available. Noticeably, the £1 = 9 estimated in this work 
is smaller than that of T N (Cl = 28) and larger than that 
of Tmi (0 = 1). This is an evidence that the pressure de- 
pendence of Tmi derived from our measurements depends 
not only on the probable multipole OP in phase II but 
also on the magnetic OP in phase III. Instead, we think 
that separating Tn and Tmi under pressure is much more 
difficult than at ambient pressure due to the enhanced 
interplay between the two phases. In addition, the origin 
of the enormous f2(T N ) may be approachable by HF be- 
haviors. A large Griineisen parameter is usually observed 
in HF systems j33| reflecting the large volume effect of 
the quasi p-f hybridization band. Indeed, Kondo behav- 
iors (e.g., the — InT dependence of electrical resistivity) 
are observed in SmRii4Pi2 below 50 K 3] and elastic 
softening at low temperatures (below 3 K) similar to HF 
behavior is also observed |23l |. 

A strong interplay between phases II and III suggested 
by the large £1 hints at an octupole scenario for Tmi in 
SmRu4Pi2. This is clear because an octupole is mag- 
netic, the same as the dipole OP in phase III, while 
the quadrupole is non-magnetic and holds TRS. In fact, 
a mix between OP of phase II (octupole) and OP of 
phase III (dipole) is proposed to explain the obscurity 
of the II-III phase transition in ref. [l4| . Noticeably, en- 
hanced elastic softening under pressure indicates that 
quadrupole moment may be one of the possible OPs in 
phase II or III under hydrostatic pressure. We infer two 



possibilities for this enhanced elastic softening. First, 
quadrupole ordering is induced by pressures in phase 
II other than the primary OP, so \gr\ and thus elastic 
softening is enhanced by quadrupole fluctuation. Sec- 
ond, quadrupole moment is an accompanying or induced 
OP in phase III besides the magnetic dipole. Therefore, 
if Tn is really approaching Tmi as anticipated from the 
Griineisen parameter analysis, the second inference will 
also produce enhanced elastic softening. There are three 
kinds of multipole ordering available in the r 6 7 ground 
state |33l ]. i.e., magnetic dipole, electric quadrupole, and 
magnetic octupole. Of these, octupole ordering may in- 
herently induce quadrupole ordering ingredients, as in 
Np02 18], or may induce quadrupole ordering by the 
mode-mixing effect |l4i l34j . or by other mechanisms. 
Thus, induced quadrupole OP in phase II or III sup- 
ports an octupole scenario for the MI transition. On 
the other hand, the largest elastic softening at around 
0.5-0.6 GPa means a kind of criticality between the two 
transitions at this pressure, considering the very differ- 
ent Griineisen parameters. In fact, the virtual tempera- 
ture dependences of Tmi and Tn based on the Griineisen 
parameters (broken lines in Fig. 0) cross at the same 
pressure. When increasing the pressure, the system pre- 
sumably moves away from the criticality state and conse- 
quently the quadrupole-strain coupling and quadrupole 
fluctuation weakens again. This argument needs to be 
checked theoretically. Here it should be noted that we 
did not observe clear signs that T N exceeds T M i under 
pressure. So the trends of the two phase transition tem- 
peratures at pressures higher than the criticality remain 
an open question. 



V. SUMMARY 

To summarize, we measured the longitudinal elastic 
constant Cl and transverse elastic constant Cr for poly- 
crystalline SmRu4Pi2 under hydrostatic pressure. Signif- 
icant elastic softening is induced by hydrostatic pressure, 
especially in the transverse Ct ■ The results give evidence 
of degeneracy of the CEF ground state with respect to 
quadrupole moments. The enhanced elastic softening in- 
dicates that quadrupole moment may be one of the pos- 
sible OPs in phase II or phase III under pressure. The 
MI transition is found to increase by applying pressure 
with a large Griineisen parameter of 9. We also found 
that elastic softening above Tmi has a strange pressure 
dependence. The greatest elastic softening is observed 
at around P = 0.5-0.6 Gpa, where the Jahn- Teller energy 
Ejt shows the largest values and |Tq| decreases. This 
strangeness may be related to the enormous Griineisen 
parameter of Tn in contrast to that of Tmi, and sug- 
gests a strong interplay between phase II and III. Our 
results prefer the octupole moment as the primary OP in 
phase II. Experiments in magnetic fields in addition to 
hydrostatic pressures are currently in progress, in order 
to separate the two successive transitions and establish a 
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detailed P-H-T phase diagram. the Ministry of Education, Culture, Sports, Science and 

Technology of Japan. 
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